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ABSTRACT 

We report the detection of a strong, organised magnetic field in the 09IV star 
[ HD 57682, using spectropolarimetric observations obtained with ESPaDOnS at the 

f*^) ■ 3.6-m Canada-France-Hawaii Telescope within the context of the Magnetism in Mas- 

sive Stars (MiMeS) Large Program. From the fitting of our spectra using NLTE model 
atmospheres we determined that HD 57682 is a 17tg 9 M star with a radius of 
Q^ \ 7.0j^g R , and a relatively low mass-loss rate of 1.4+gg 5 x 10~ 9 M Q yr -1 . The pho- 

. tospheric absorption lines are narrow, and we use the Fourier transform technique to 

infer v sin i — 15 ± 3 kms -1 . This v sinz implies a maximum rotational period of 31.5 
d, a value qualitatively consistent with the observed variability of the optical absorp- 
■ tion and emission lines, as well as the Stokes V profiles and longitudinal field. Using 

?—i ' a Bayesian analysis of the velocity-resolved Stokes V profiles to infer the magnetic 

field characteristics, we tentatively derive a dipole field strength of 1680l3gg G. The 
derived field strength and wind characteristics imply a wind that is strongly confined 
by the magnetic field. 

Key words: stars: magnetic fields - stars: winds - stars: rotation - stars: early-type 
- stars: individual: HD 57682 - techniques: spectropolarimetry 



1 INTRODUCTION 

The phenomenon of magnetism in hot, massive OB stars 
is not well studied. To date, repeated detections of circular 
polarization within line profiles have firmly established the 
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presence of magnetic fields in three O-type stars (6 1 Ori C, 
HD 191612, and ( Ori A; Donati et al. 2002, 2006; Bouret 
et al. 2008) and a handful of early B-type stars (e.g. Do- 
nati et al. 2001; Neiner et al. 2003; Alecian et al. 2008; Petit 
et al. 2008; Silvester et al. 2009). Indications of magnetic 
fields have also been reported in a number of other OB stars 
(e.g. by Hubrig et al. 2008; 2009). Such objects represent vi- 
tal targets for the study of stellar magnetism. Their strong, 
radiatively-driven winds couple to magnetic fields, gener- 
ating complex and dynamic magnetospheric structures (e.g. 
ud-Doula & Owocki 2002), which modify mass loss, and may 
enhance the shedding of angular momentum via magnetic 
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braking (e.g. ud-Doula et al. 2008). The presence of even 
a relatively weak magnetic field can profoundly influence 
the evolution of massive stars and their feedback effects, 
such as mechanical energy deposition in the ISM and super- 
nova explosions (e.g. Ekstrom et al. 2008). For these rea- 
sons, the Magnetism in Massive Stars (MiMeS) Large Pro- 
gram (Wade et al. 2009; Grunhut et al. 2009) is exploiting 
the unique spectropolarimetric characteristics of ESPaDOnS 
at the Canada- France-Hawaii Telescope (CFHT) and NAR- 
VAL at the Telescope Bernard Lyot (TBL) to obtain criti- 
cal missing information about the poorly-studied magnetic 
properties of these important stars. 

One of the MiMeS targets is HD 57682, an 09 sub-giant 
star. HD 57682 is reported to be a runaway star (Comeron et 
al. 1998) with no known companions (e.g. de Wit et al. 2005; 
Turner et al. 2008) and shows no photometric variability 
(e.g. Balona 1992). However, HD 57682 exhibits variability 
of UV lines characteristic of magnetic OB stars (e.g. Schnerr 
et al. 2008), and this is the primary reason it was included 
in the MiMeS target list. 



2 OBSERVATIONS 

Eleven high- resolution (R ~ 65,000), broadband (370- 
1050 nm) circular polarization (Stokes V) spectra of 
HD 57682 were obtained with the ESPaDOnS spectropo- 
larimeter, mounted on the 3.6-m CFHT, as part of the Sur- 
vey Component of the MiMeS Large Program. The spectra 
were acquired and reduced in a manner essentially identical 
to that described by Silvester et al. (2009). The observations 
were obtained during two separate ESPaDOnS runs. The 
first two observations were obtained on successive nights in 
2008 Dec. (high extinction hindered our first observation), 
and an additional 9 observations were obtained in 2009 May 
on five different nights. On nights when more than one ob- 
servation was obtained, the un-normalised spectra were co- 
added prior to analysis, yielding 7 independent observations; 
A 5 th -order polynomial (or lower) was used for normalisa- 
tion of each individual order. A complete summary of the 
ESPaDOnS observations is given in Table [T] 

Least-Squares Deconvolution (LSD; Donati et al. 1997) 
was applied to all polarimetric observations in order to in- 
crease the signal-to-noise ratio (S/N) for this analysis. We 
adopted a mask in which all the Stark-broadened hydrogen 
and helium lines, and the metallic lines that are blended 
with the hydrogen lines were excluded. This resulted in ap- 
proximately 165 relatively pure photospheric lines used for 
the creation of the mean, normalised Stokes /, V, and diag- 
nostic null (N) profiles, yielding an increase in the Stokes V 
S/N by a factor of ~7.5. All LSD profiles were computed 
on a spectral grid with a velocity bin of 3.6 kms -1 . The 
Stokes V, diagnostic null, and the Stokes / profile of our 
2008 Dec. 6 observation are shown in Fig. \T\ (a compilation 
of all Stokes / and V profiles is provided in Fig. [5}. A careful 
examination of the corresponding reduced spectrum clearly 
shows the presence of Stokes V signatures of similar mor- 
phology in numerous individual absorption lines of different 
chemical species. We also list the false alarm probabilities 
according to the criterion of Donati et al. (2002; 2006) in Ta- 
bleQ] resulting in 5 definite detections, 1 marginal detection, 
and 1 non-detection. 
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Figure 1. LSD Stokes V (top), null N (middle), and Stokes / 
profiles of HD 57682 from 2008 Dec. 6. The V and N profiles arc 
expanded by the indicated factor and shifted upwards for display 
purposes. A clear Zeeman signature is detected in the Stokes V 
profile, while the null profile shows no signal. The integration 
limits used to measure the longitudinal field are indicated by the 
dotted lines. 

The longitudinal magnetic field (Be) was inferred from 
each LSD Stokes V profile in the manner described by Sil- 
vester et al. (2009). The longitudinal field measurements, 
which are reported in Table [T] vary between -167 G and 
+282 G, with a typical uncertainty of 20-40 G. While the 
longitudinal field provides a useful statistical measure of 
the line-of-sight component of the field, we stress that we 
do not use it as the primary diagnostic of the presence of 
a magnetic field. This is because a large variety of mag- 
netic configurations can produce a null longitudinal field. 
However, nearly all of these configurations will generate a 
detectable Stokes V signature in the velocity-resolved line 
profiles, which is what we observe. 

In addition to the ESPaDOnS observations, archival 
IUE data were used to constrain the fundamental param- 
eters of HD 57682 (see Sect. and investigate its UV and 
optical line variability (see Sect. [3]). Three IUE spectra were 
obtained over a period of 5 years, one on 1978 Dec. 12 and 
the other two on 1983 Nov. 20. The first observation was 
obtained with the small aperture, while the latter were with 
the large aperture. 



3 STELLAR PARAMETERS 

Non-LTE, expanding atmosphere models calculated with the 
CMFGEN code were used in our analysis (see Hiller & Miller 
1998 for details) . A summary of the fundamental parameters 
obtained is reported in Table [2] The surface gravity (logg) 
and effective temperature (T e ff) were derived from the opti- 
cal spectrum. The H<5 and H7 wings provided a measurement 
of log g, while the He I 447- nm and He 11 454-nm transitions 
were used as the main diagnostic for T e g. However, fits to 
other helium transitions were also checked for consistency. 
Model fits to the He 11 454-nm and H7 are shown in Fig. [2] 
The uncertainties quoted in Table [5] already take into ac- 
count the observed temporal variability (see Sect. [4}- The 
stellar luminosity L* was assumed to be typical of other 
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Table 1. Journal of ESPaDOnS observations listing the date, 
the heliocentric Julian date (2,454,000+), the number of sub- 
exposures and the exposure time per individual sub-exposure, 
and the peak signal-to-noise (S/N) ratio (per 1.8 kms" 1 velocity 
bin) in the Null spectrum in the V-band, for each night of obser- 
vation. Columns 5-6 list the false alarm probability and the mean 
longitudinal field inferred from the LSD profiles. 



Table 2. Summary of stellar and wind properties of HD 57682. 



Date 


HJD 


£cxp 


S/N 


FAP 




B e ± o B 






(fl) 








(G) 


2008/05/12 


806.080 


4 x 500 


259 


5 x 10~ 


4 


282 ±114 


2008/06/12 


806.786 


4 x 500 


937 


< 10" 


8 


266 ± 30 


2009/04/05 


955.768 


8 x 600 


1502 


3 x 10" 


7 


-46 ± 21 


2009/05/05 


956.750 


8 x 540 


1401 


1 x 10" 


7 


-35 ± 22 


2009/07/05 


959.102 


4 x 540 


667 


6 x 10" 


2 


-94 ± 50 


2009/08/05 


959.749 


8 x 540 


1136 


1 x 10- 


8 


-167 ±30 


2009/09/05 


960.748 


8 x 540 


915 


6 x 10- 


7 


-109 ±37 
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Figure 2. Modelling of the He II 454-nm line with atmospheric 
models corresponding to different T c ff values (left panel; log g = 
4.0), and the H7 line with models corresponding to different logg 
values (right panel; T off = 34, 500 K) for the night of 2009 May 
4. Note that only the wings of the H7 profile are used; the core 
is likely contaminated by circumstellar material (see Sect.|6]l. 



Spectral type 

T cff (K) 

log 9 (cgs) 

R* (R Q ) 

log(i*/L Q ) 

M* (M Q ) 

log M (Moyr- 1 ) 

Voo (kms" 1 ) 

l°g(£x/T Soi ) 



09IV (Walborn 1972) 

34 500 ± 1000 

4.0 ± 0.2 

7 + 2 - 4 
'■ u -1.8 

4.79 ± 0.25 
17 + 19 

11 -9 

-8.85 ± 0.50 

1200 + 500 
izuu _ 200 

-6.34 




Velocity (km/s) 

Figure 3. Temporal variations of the Ho (left) and He I 447-nm 
(right) lines throughout our observing run. The time-averaged 
profile is plotted in red (dashed) to emphasise variations. The 
date of observation is indicated in the Ha frame. 



Due to the lack of wind diagnostic lines, clumping was not 
included in our analysis. Auger ionizations from X-rays were 
taken into account in the models, since they are essential in 
stars with spectral types similar to HD 57682 (see Martins et 
al. 2005b). The level of X-ray emission used in our CMFGEN 
models is taken from ROSAT measurements (Berghoefer et 
al. 1996), scaled to our distance. 



late O-type stars (see Martins et al. 2005a, Table 4) . There- 
fore, through the analysis of IUE spectra we derived a dis- 
tance of ~ 1.3 kpc and a corresponding reddening parameter 
E(B - V) of 0.07 from the fit to the observed UV contin- 
uum. The obtained values are in good agreement with other 
modern literature estimates (e.g. Wegner 2002). The stellar 
radius was directly computed from R* = (L^/AnaT^^) 1 ^ 2 
and the stellar mass follows from = gR^/G. Standard 
errors for R t and M* were computed assuming an error in 
the luminosity of ±0.25 dex (see Martins et al. 2005b). Test 
models were also investigated with different abundances for 
the CNO elements, but precise values were difficult to de- 
termine and we found no clear evidence of chemical enrich- 
ment or depletion with respect to solar abundances of either 
Grevesse & Sauval (1998) or the more recent Grevesse et al. 
(2007). A deeper study is necessary, but beyond the scope 
of this paper. 

The wind parameters M and Voo were determined from 
the fit to the C IV 154.8, 155.1-nm profile, from IUE data 
SWP21590, and appears consistent with all IUE data avail- 
able. A standard /3-velocity law was considered, with (3 = 1. 



4 TEMPORAL VARIABILITY AND 
ROTATION 

Our first observations obtained in 2008 Dec. show an en- 
hanced Ha emission and stronger absorption line depth in 
other elements compared to the observations obtained in 
2009 May (see Fig. [3] for an example with He lines) . The 
stability of neighbouring lines and the use of a low-order 
polynomial during normalisation imply that the variability 
is real. An analogous anti-correlation is seen in the magnetic 
O-type star 8 1 Ori C (Donati et al. 2002; Wade et al. 2006). 
We also note that the longitudinal magnetic fields measured 
in 2008 Dec. versus 2009 May show opposite signs (see Ta- 
ble [TJ . Additionally, the IUE U V spectroscopy shows sig- 
nificant variations of the 139.4, 140.2-nm Si IV doublet (see 
Fig. [4j - a phenomenon that is not observed in other UV 
lines and occurs exclusively in magnetic B-type stars, e.g. 
£ Cas (Neiner et al. 2003) and all magnetic He-strong stars. 
Our CMFGEN models are not able to reproduce this vari- 
ability by increasing the mass-loss rate, which, as we show 
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1386 1389 1392 1395 1398 1401 1404 1407 




Velocity (km s ) 

Figure 4. Overplot of the Si IV line profiles of the three high- 
resolution IUE spectra of HD 57682. The two doublet rest wave- 
lengths are indicated by vertical dashed lines. The lower part dis- 
plays the significance of the variability, expressed as the square 
root of the ratio of the measured to the expected variances. Note 
that the variability observed at 1398 A is well-known and due to 
small differences in the echelle order overlap corrections. 

in Sect. [6] leads us to believe that it is due to magnetically- 
confined circumstellar material. 

In order to determine the rotational velocity usini of 
HD 57682 we computed the Fourier transform of many un- 
blended absorption lines throughout the optical spectrum 
(e.g. Gray 1981). The results of our analysis suggest that 
there is an important non-rotational contribution to the 
line profile. Our CMFGEN models employed an additional 
40 kms -1 macro-turbulence to provide a better fit to the 
observed absorption lines. From the analysis of about 30 
different lines from the 2008 Dec. 6 and 2009 May 4 obser- 
vations, we infer a mean vs'mi of about 15 ± 3 krns" 1 , in 
good agreement with Balona (1992). 

Due to the limited temporal sampling of our observa- 
tions, we are unable to unambiguously determine the rota- 
tional period of HD 57682. However, using an upper limit 
on our inferred radius of 9.4 Rq and the inferred rotational 
broadening of 15 kms" 1 , we estimate that the maximum 
(rigid) rotation period P rot is equal to about 31.5 d. The ob- 
served variations of Ha (see Fig. [3} are qualitatively consis- 
tent with such a period, because the night-to-night changes 
are small compared to the difference between the mean pro- 
file corresponding to the two different observing runs. As 
well, the variability of the Stokes V profiles and longitudi- 
nal field measurements also support a variability timescale 
considerably longer than 5 nights. We therefore conclude 
that the variability of the magnetic and spectroscopic ob- 
servables are consistent with rotational modulation, with a 
period of a few weeks. 



5 MAGNETIC FIELD 

The longitudinal field measurements imply a strength of the 
dipole component of the magnetic field of at least 1 kG, con- 
sistent with the strengths of fossil fields studied in Ap/Bp 
stars, as well as other B and O-type stars. To estimate the 
surface magnetic field characteristics of HD 57682, we used 
the same method as Petit et al. (2008; and in prep.), which 
compares the observed Stokes V profiles to a large grid of 
synthetic profiles, described by an oblique-rotator model, to 
estimate the magnetic properties of the star. The model is 
parametrized by the dipole field strength B p , the rotation 



axis inclination i, the positive magnetic axis obliquity /3, 
and the rotation phase <p. In order to compute synthetic 
Stokes V profiles, we need to infer the line properties using 
the Stokes / profile. As the detailed formation of the line 
profile and the source of its variation is uncertain, we opted 
for the simplest assumption and used the mean of the ob- 
served profiles to determine the intensity profile parameters. 

Assuming that only ip may change between different 
observations of a given star, the goodness-of-fit of a given 
rotation-independent (B p , i, /3) magnetic configuration can 
be computed to determine a configuration that provides the 
overall best likelihood for all the observed Stokes V profiles. 
In Fig. [5] we compare the synthetic profiles of the overall 
best likelihood configuration (in green) with the best likeli- 
hood configuration for each observation (in red). The quality 
of the different fits are similar, showing that a single inclined 
dipole can reproduce the observations as well as individual 
dipole configurations, although both are imperfect. However, 
any features that cannot be explained are treated formally 
as additional noise. 

Fig. [5] shows the resulting marginalised a posteriori 
probability density function (PDF) for the dipole field 
strength (bottom left). The inferred value for the magnetic 
surface field strength in the 68.3% credible region, is then 
1680 

-356 The value of the dipole field strength will vary 
greatly with a small change of the inclination value, as il- 
lustrated in the B p -i 2D PDF shown in Fig. [5] (top left). A 
determination of the rotation period by further observations 
will improve the estimate of the dipole field strength. 

Based on these data, the obliquity of the magnetic axis 
is closely correlated with the inclination as well, as demon- 
strated by the (3-i 2D PDF shown in Fig. [5] (top right). 
According to the 99.7% credible regions of the marginalised 
PDFs of i and f3 (shown in Fig. [5] in blue and red respec- 
tively), we infer that the obliquity is between 10° and 50° if 
the inclination is between 47° and 84° , or between 130° and 
169° if the inclination is between 95° and 132° . 



6 DISCUSSION AND CONCLUSIONS 

It is interesting to note that the mass-loss rate derived from 
the Ha profiles is inconsistent with mass-loss derived from 
the C IV UV lines. Test models indicate that the intensity of 
the Ha emission observed at different dates requires mass- 
loss rates above ~10~ 7 Mq yr _1 . With such high values, the 
predicted UV spectra would present strong P-Cygni profiles 
(e.g. Marcolino et al. 2009), which are not seen in any of 
the IUE spectra available for HD 57682. In fact, HD 57682 
presents a UV spectrum typical of Galactic weak wind stars, 
with no significant P-Cygni emissions (Martins et al. 2005b; 
Marcolino et al. 2009). This leads us to believe that the Ha 
profiles are not formed in the wind, but are possibly a result 
of confined circumstellar material. Walborn (1980) already 
reported this line to be unusual and of non-nebular origin. 

A preliminary Bayesian analysis suggests that the mag- 
netic field of HD 57682 is approximately reproduced by a 
dipole field with a characteristic surface field strength of 
1680 Iggg G. We computed the magnetic wind confine- 
ment parameter (ud-Doula & Owocki 2002) 1.4 x 10 4 for 
HD 57682, using B p ~ 1600 G, and the physical param- 
eters listed in Table [2] (77, ranges from 4 x 10 3 - 2 x 10 4 
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Figure 5. Left: Marginalised a posteriori probability density for the dipole strength (bottom left), the magnetic field obliquity (bottom 
right and solid-filled red, with the inclination PDF overplotted in line-filled blue, the B p -i plane (top left) and the i-f3 plane (top right). 
The 68.3%, 95.4% and 99.7% credible regions are indicated in dark to light shades respectively. The additional contours for the 2D planes 
are 99.994% and 99.9999%. Right: Our time series in Stokes V are shown, along with individual geometries that provide the maximum 
likelihoods for single observations (in thick green) and with the geometry that provides the maximum likelihood for all the observation 
combined (in thin red). The associated reduced x 2 are a l so indicated. 



using the 99.7% credible regions). Although strongly con- 
fined, this places HD 57682 at an interesting intermedi- 
ate regime amongst OB stars, between the super-strongly- 
confined wind of a Ori E {rj„ ~ 10 7 ) and the much more 
weakly confined wind of 6 l Ori C (ry» ~ 20). Magnetic con- 
finement of the wind of HD 57682 would naturally explain 
the Ha variability, the UV variability, and the line depth 
variability. Moreover, shocks produced in a magnetically 
channelled wind could boost the X-ray emission, providing 
an explanation of the high X-ray luminosity without requir- 
ing a high mass- loss rate (e.g. Cohen et al. 2008). However, 
it is puzzling that HD 57682 presented a soft X-ray spectrum 
during the ROSAT observations, contrary to what would be 
expected from a fast wind and the strong magnetic confine- 
ment. 

This discovery brings the number of firmly-established 
magnetic O-type stars to four. With its strong wind confine- 
ment, and its stellar and wind properties determined, HD 
57682 provides a testbed for future MHD simulations. How- 
ever, more data are needed to better determine the surface 
dipole field strength, magnetic obliquity, and inclination. 
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